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I. INTRODUCTION 



The investigation of heavy flavour (charm or bottom) production in proton-proton col- 
hsions at the LHC (Large Hadron Colhder) are important for testing perturbative QCD 
calculations in a new energy domain, where very small Bjorken-x momentum fractions are 
expected to be probed. The detection of the heavy hadrons containing charm or bottom 
quarks can be done in several ways, either through their non-leptonic weak decays or through 
their semileptonical decays, where electrons or muons are measured inclusively. 

Already at the RHIC (Relativistic Heavy Ion Colhder), the PHENIX and STAR Collab- 
oration |2| measured the production of muons and electrons from heavy flavour decays in 
pp collisions at = 0.2 TeV. These data have been compared with perturbative QCD cal- 
culations in the FONLL framework jsl and were found in agreement with the measurements 
within the experimental and theoretical uncertainties. 

Experimentally the measurement of the semileptonic cross sections suffers from a large 
background due to the semileptonic decay of primary light hadrons including pions and 
kaons (the main contribution) and other meson and baryon decays (such as J/ip and low 
mass resonances rj, p, u and 0), from secondary leptons produced from secondary light 
hadron decays, and from secondary hadrons escaping from material surrounding the tracking 
chambers. All these backgrounds must be subtracted based on Monte Carlo simulations 
using the usual event generators to obtain the semileptonic yield due to the decay of heavy 
hadrons. 

The measured cross sections for the inclusive production of leptons from heavy hadron 
decays are of two types. In two experiments the production cross section is measured 
separately for solely 6-hadron decays ^, |5|. In the CMS experiment [4| the 6-hadron cross 
section was discriminated by measuring the muon transverse momentum with respect to 

n 

the closest jet. In the ALICE experiment |5[ the production cross section of electrons, 
(e"*" + e~)/2 from semileptonic bottom hadron decays was selected by using the information 
on the distance of the secondary decay vertex with respect to the primary vertex. Due to 
their long lifetime bottom hadrons decay at a secondary vertex displaced in space from the 
primary collision vertex. In all other measurements of the semileptonic production which 
were presented by the ATLAS [6] and ALICE jr-Q] collaborations no attempt was made to 
separate leptons from charm and bottom hadrons. Therefore these cross sections constitute 
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the sum of lepton (electrons or muons) production cross sections of all charmed hadrons 
{D^, D^, Ds and c-baryons) and all bottom hadrons {B^, B^, Bg and 6-baryons). 

At LHC energies the cross sections of charm and bottom production using special non- 



eptonic decays of charm and bottom hadrons have been reported. So, ALICE 10|], ATLAS 
ll| and LHCb ^Ji] presented cross section data on , D^, D*^ and Dg production and 
CMS [l^ published measurements for B^, Bg and A;, production cross sections. All these 
cross sections have been calculated in the framework of the general- mass- variable-flavour- 
number-scheme (GM-VFNS) [l^, which has been developed in the last ten years for various 
processes. In this scheme we calculated the inclusive B meson production cross section in 
pp collision [15] at = 1.96 TeV and in pp collisions 16] and found good agreement 
with the respective data from the CDF run II 171] and also with the data from the CMS 



collaboration 



13 1 at the LHC at = 7 TeV. Besides the data based on special non-leptonic 



decays of B mesons mentioned above there exist also several cross section measurements for 
pp — )■ BX and pp — )■ BX followed by the inclusive decays B — )• J /ipX and B — )• il){2S)X. 
These cross sections have been calculated recently also in the GM-VFNS 18] a.nd compared 
with experimental data from CDF 17, 19], CMS Q, LHCb 2li, ATLAS 



m and ALICE 



231 ] collaborations at the LHC. The agreement between the calculated cross sections and 



the data was satisfactory. In the same scheme the inclusive production of charmed mesons 
D^,D~^,D*~^,Ds and the charmed baryon A,, has been calculated and the results for the 
ALICE |l^, ATLAS H] and LHCb 12] kinematic conditions have been presented in ref. 
24] and compared to the experimental data in the respective presentations of the ALICE 
10], ATLAS [nl and LHCb 12] collaborations and for the ALICE data also in ^ . 

The GM-VFNS is similar to the zero- mass variable- flavour- number scheme (ZM-VFNS), 
in which the heavy quark m is neglected in the calculation of the hard-scattering cross 
sections. The predictions in the ZM-VFNS are expected to be reliable only in the region of 
very large values of the transverse momentum pt of the produced heavy hadron since terms 
of the order of m? /p^ are neglected in the hard-scattering cross sections. In the GM-VFNS 
these /p^ terms are retained as they appear in the so-called fixed-fiavour-number scheme 
(FFNS) in such a way that by applying appropriate subtractions to the FFNS the GM- 



VFNS approaches the ZM-VFNS with the usual MS prescription in the limit pt/ttl oo. 
Whereas in the FFNS the gluon and the light partons are the only active partons in the initial 
state and the heavy quark appears only in the final state, produced in the hard-scattering 
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process of light partons, in the GM-VFNS as in the ZM-VFNS, due to the subtraction of 
the mass-singular contributions in the initial and final state, the heavy quark appears also 
in the initial state, and for the final state appropriate fragmentation functions (FFs) for the 
transition heavy — quark — )■ heavy — hadron must be introduced which absorb the collinear 
singular contributions in the final state. Such FFs for the transitions b — > -B, b,c ^ D^, 
and D*^ and b,c ^ and have been extracted in ref. 



15 



25 



26|] at NLO in the 



MS factorization scheme with nj = 5 (for b) and with np = A (for c) flavours, consistently 
within the GM-VFN framework using data for the scaled energy (x) distribution da/dx of 
e+e" BX and e+e" DX etc, respectively, measured by the CERN LEPl, SLAG SLG, 
GESR GLEG and the KEKB Belle collaborations. 

The content of this paper is as follows. In Sec. 2 we summarize our input choices of 
PDFs and B- and D-meson FFs. In this section we also explain how the fragmentation of 
these mesons into leptons has been calculated. In Sec. 3 we compare the predictions of the 
GM-VFN scheme with the existing data from the recent LHG run at \/S = 7 TeV 4i9\ and 
= 2.76 TeV Qj. We end with a summary in Sec. 4. 



II. SETUP, INPUT PDFS AND FFS 



The theoretical framework and results of the GM-V^ 



have been previously represented in detail in refs. 14 



'N 



approach for pp {pp) collisions 



15| . Here, we describe our choice 



of input for the numerical analysis for inclusive lepton (e or yu) production from charm and 
bottom hadrons. We use for the ingoing protons the PDF set GTEQ 6.6 as implemented 
in the LHAPDF 28|| library. This PDF set was obtained in the framework of a general-mass 
scheme using the input mass values rric = 1.3 GeV, and rrify = 4.5 GeV, and for the QGD 
coupling ai^\mz) = 0.118. The c-and 6-quark PDFs have the starting scales /xq = "^c and 
/io = rrib, respectively. 

The nonperturbative FFs for the transition of b — t- B were obtained by a fit to e^e" 
annihilation data from the ALEPH 



been presented in 



29|, OPAL [301] and SLD [3l| collaborations and have 



15| . The combined fit to the three data sets was performed using the 

0.1181 adopted from 



NLO _scale parameter A^^ = 227 MeV corresponding to af\mz 



ref. [27|. Gonsistent with the chosen PDFs, the starting scale of the b ^ B FF was assumed 
to be /io = TTLh, while the q,g ^ B FFs, where q denotes the light quarks including the charm 
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quark, were taken to vanish at /io- As input we used the FFs with a simple power ansatz 
which gave the best fit to the experimental data. The data from OPAL and SLD included all 
6-hadron final states, i.e. all B mesons, B^, B^ and B^ and 6-baryons while in the ALEPH 
analysis only final states with identified B^ {B^) and B^ (5°) were taken into account. For 



the fit in ref. 



15| it was assumed that the FFs of all 6-hadrons have the same shape. This 



will also be assumed in this calculation. In addition we assumed that all 6-hadrons have the 
same branching fractions and decay distributions into leptons as one of the B mesons, B^ or 
B^ . The only difference results from the different fractions b b — hadron, which are taken 
from the Particle Data Group [32]. Based on these values the prediction, for example for 
B^, is multiplied by 2.49. The bottom mass in the hard scattering cross sections is rUb = 4.5 
GeV as it is used in the PDF CTEQ6.6 and in the FFs for b B etc. 

For the transitions c — )■ D^, we employ the FFs determined in 25|. They are based on 
: its to the most precise data on D meson production from the CLEO Collaboration at CESR 
33| and from the Belle Collaboration at KEKB 3^. Actually there are several alternative 
fits presented in 25|. Here we use the so-called Global-GM fit, which includes fitting in 
addition to OPAL data 35 1 together with the CLEO and Belle results. The fits in I25I are 
based on the charm mass rric = 1.5 GeV, which is slightly larger than the one used in the 
CTEQ6.6 PDF fits. The starting scale for c — )■ -D is /iq = rric as it is for the g,q D FFs, 
whereas for the b ^ D FF it is /io = iTib- 

The subtractions related to renormalization and to the factorization of initial- and final- 
state singularities require the introduction of scale parameters fiR, /i/ and fip- We choose 
the scales to be of order m-r, where rriT is the transverse mass rriT = \f^^^rrr? and m = rrib 
for the case of bottom quark production and m = rric for charm quark production. For 
exploiting the freedom in the choice of scales we introduce the scale parameters C,i (i=R, I, 
F) by setting /ij = C,i ttlt- To describe the theoretical uncertainties we vary the values of 
independently by factors of two up and down while keeping any ratio of the parameters 
smaller than or equal to two. The uncertainties due to scale variation are the dominating 
source of theoretical uncertainties. Therefore PDF related uncertainties and variations of 
the bottom and charm mass are not considered. 

The fragmentation of the final state partons i into lepton I {I = or /x^) is calculated 
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from the convolution 



Jx z z Tb dz 



(1) 



In this formula Di^B{x, fJ-p) is the nonperturbative FF determined in 15| for the transition 



i ^ B and in 



25| for i D (with subscript B replaced by in (1)), is the total B decay 



width {Tb — ^ Fd in case oi i ^ D) and dV{z, Pb) I dz is the decay spectrum of 5 — i- / or 
D ^ I, respectively. For given lepton transverse momentum pt and rapidity y, Pb is given 
by Pb = |Pb| = Pta/1 + sinK^y/z. The decay distribution dV/dk'^, where the momentum 
fc^ is parallel to Pb, is obtained from the decay distribution in the rest system of the B 
meson using the formula (3.16) in ref. 36|, where this formula was derived for the decay 
B — !• J/ip instead of 5 — )■ /. This leads to dT{z, PB)/dz used in eq.(l) with z = k'j^/PB- 

The electron energy spectrum in inclusive B — t- evX decays has been measured by the 
BABAR Collaboration |37| in the range > 0.6 GeV up to its kinematical limit. This 
spectrum has been measured at the T(45') resonance. The partial branching fraction for 

> 0.6 GeV has been determined as B[B evX^E^ > 0.6 GeV] = [10.36 ± 0.06(stat) ± 
0.23{syst)]% 

Following an ansatz for the momentum spectrum as a function of p in GeV as given in 



ref. 



we fitted the BABAR spectrum in the B rest system using the formula (3.9) in 
361] for the transformation of the spectrum in the B rest system to the BABAR laboratory 
frame, where |Pb| = 0.341 GeV by the following formula 

-24.23 + 14711.2 exp[-l. 73 ln2[2.74 - l.lOp]]) 



(2) 



(-41.79 + 42.78 exp[-0.5( ^(p - 1.27)71.8)^-^^]) 
This differs from the fit in 



38| which was done in the BABAR T(4S') laboratory frame. 



/b(p) is related to the partial semileptonic decay spectrum dT/dp by dV /dp = cfB^p) with 



c = 4.27151 10^^ according to j36|, where dT/dp has units [ps'^GeV^^]. The quality of the fit 
can be seen in Fig.l (left figure) together with the effect of the boost into the B rest system, 
which indeed is small. This fit yields for the partial branching fractions B{B — )■ evX^p > 
0.6 GeV) = 11.04% and B{B — )■ evX,p > 0) = 12.35%. The latter branching fraction agrees 
approximately with the more recent PDG value B{B — t- evX^p > 0) = [10.74 ± 0.16]% |32|. 



The inclusive electron spectra for the decays — )■ e^vX and — )■ e^vX have been 



measured by the CLEO Collaboration 



on the ?/'(3770) resonance (75 = 3.73 GeV) 




p[GeV] P[GeV] 

FIG. 1: fit for inclusive lepton spectrum of B decays to BABAR (left frame) data in te BABAR 
laboratory frame together with spectum boosted to the B rest frame and to inclusive lepton 
spectrum of D decays to CLEO (right frame) data in the CLEO laboratoty frame together with 
spectrum boosted to the D rest frame. 

as a function of the electron momentum p. The measured partial branching fractions for 
p > 0.2 GeV are B{D+ e+uX,p > 0.2 GeV) = [14.97 ± 0.19{stat) ± 0.27{syst)]% 
and B{D° euX,p > 0.2GeV) = [5.97 ± 0.15(stat) ± 0.10{syst)]% from which, on the 
basis of fits to both spectra by the CLEO Collaboration, the following branching fractions 
B{D+ e+iyX,p > 0) = [16.13 ± 0.20{stat) ± 0.33{syst)]% and B{D° e+uX,p > 0) = 
[6.46 ± 0.17{stat) ± 0.13{syst)]% result. The corresponding PDG values are: B{D^ — >■ 
e+uX) = [16.07 ± 0.30]% and B{D^ e+uX) = [6.49 ± 0.11]% H, which agree vey well 
with the CLEO values. The spectra of and D'^ as measured by CLEO coincide very well 
inside errors in the measured p range 0.2 < p < 1 GeV. In the CLEO laboratory system the 
momentum |Pd| is 0.243 GeV and the momentum is 0.277 GeV. These numbers are 
needed for the fits in the and rest systems. The result of the fit for the spectrum in 
the D rest system, for example for the it is 

/^(p) = 17.91(p + 0.0034)2 *56(0.98 - p)^-^^ (3) 

fnip) in eq.(3) is equal to the lepton spectrum foip) = dV /dp in units of [ps^^GeV^^]. The 
quality of the fit to the CLEO data is shown in Fig{T] in the right frame together with the 
spectrum in the D rest system. As can be seen there the difference between the two spectra 
is really small. The branching ratio for this spectrum is: B{D^ — t- e'^ + uX) = 15.29%, which 
is sufficiently close to the valuues from CLEO given above. With these parametrizations 
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of the electron spectra in the B, respectively D, rest system, we calculated the lepton 
spectra dT/dx in the moving system as function of /c^ = xPn (xPd) where k'j^ is the lepton 
momentum that is parallel to Pb(Pd) using eq.(3.17) in 36| . 

To simplify the calculation we applied for dT/dx the asymptotic formula as given in 36| . 
which is the approximation for ^ Mb {Pd ^ ^^d)-, where Mb ans M^i are the masses 
of the B and D meson, respectively. We calculated dT/dx for various Pb and found that 
the exact formula differs from the asymptotic formula by less than 5% for Pb = 10 GeV. 
This is easily achieved even for the smaller pt values of the leptons in our applications 
since dT/dx is peaked at very small x, so that the average x is below 0.2 in the case of B 
decays and similarly for D decays. To make the calculation easier we calculated the lepton 
spectra for just one flavour state, the D^, and included the other D mesons, D^, Dg and the 
contributions of charmed baryons, as for example Ac, by an appropriate normalization factor. 
For this normalization factor is 2.294 and was calculated using a compilation of charm 
hadron production fractions by Lohrmann 391] (Table 3 of this reference) and the lepton 
branding fractions for D°,D+,Ds and equal to 0.0646%, 0.161%, 0.065% and 0.045%, 
respectively |32j. In case we would have based this normalization factor on production, 
the normalization factor would be slightly different. This difference is compensated by the 
different normalization of c — > and the different lepton branching ratio, which result in 
almost equal cross sections da / dpx for the production of lepton from charmed hadrons via 

or D+. 



III. RESULTS AND COMPARISON WITH LHC DATA 

In this section we collect our results for the cross sections da/dpj- as a function of pt- 
for the various LHC experiments and compare them with the published data. We start 
with the results on inclusive b-hadron production with muons at y/S = 7TeV by the CMS 
Collaboration 4| and the ALICE collaboration jsl. The CMS data are for the production of 
the sum of fi~^ + /i~ in the range from 6 < < 30 GeV, and are integrated over the rapidity 
y in the region —2.1 < y < 2.1. The comparison between the GM-VFNS predictions and the 
experimental data is shown in FigIS The data lie slightly above the full default curve (default 
scale choice) and agree better with the prediction with the maximal scale choice. The ALICE 
cross section da/dpr jsj for the production of electrons and positrons ((e^ + e~)/2) from 
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?1G. 2: GM-VFNS predictions for inclusive lepton production compared to CMS 4j] and ALICE 
^ data. 

bottom hadron decays at = 7 TeV in the rapidity range —0.8 < y < 0.8 is measured for 
smaller p-r < 8 GeV (note the different definition of the cross section of the ALICE data as 
compared to the CMS data). Our prediction in comparison with the ALICE data, is shown 
also in Figj2l where we included the data only for pt > 1.5 GeV. The agreement between 
the data and the prediction is quite good. The data with their errors lie quite well inside the 
theoretical uncertainty band, although nearer to the upper prediction from the respective 
scale choice (dashed line) at the larger px values. We remark that in the two predictions in 
Fig|2] we included also the contributions from 6 — > D — t- lepton. The FFs needed for this 
contribution were taken from 25| (global fit). These contributions are very small, of the 
order of 1 to 2 % compared to c -D contribution, depending on px, and are not visible in 
the logarithmic plots in FigJH The b D contributions are also included in all the following 
predictions. 

Next we present our predictions for the ATLAS and ALICE kinematical conditions for 
the combined cross sections da/dpT for electron (or muon) production resulting from the 
sum of charmed and bottom hadron production for = 7 TeV and one results from 
ALICE for \/S = 2.76 TeV. This latter measurement was done by the ALICE collaboration 
in order to have the reference pp cross section for their PbPb measurements at the same 
energy. In general the contributions from D {B) decays dominate at small (large) px and 
are approximately equal at ~ 4 GeV. The ATLAS data (sf consist of three different 
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FIG. 3: GM-VFNS predictions for inclusive leptons from bottom hadron production compared to 
ATLAS B data. 

measurements of da/dpT for the production of electrons (e^ + e^) as a function of px for 
7 < pr < 26 GeV integrated over the rapidity range —2.0 < y < 2.0 excluding the rapidity 
interval 1.37 < \y\ < 1.52, (ii) the production of muons (/i^ + in the same px and y 
range as in (i), and (iii) the production of muons {^'^ + ) in the pt range 4: < px < 100 
GeV and the y range —2.5 < y < 2.5. The result of our predictions is shown in FigJS] The 
agreement with the ATLAS data 6| is very good. The data points lie in all three plots with 
their errors inside the band obtained for the variation of the scale choices (upper and lower 
dashed curve). In the special plot for muon production (second ATLAS plot) the data lie 
very near to the default choice prediction (full curve) inside their small experimental errors. 
The same is true for the third plot in FigJSl Here the data and the default prediction agree 
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FIG. 4: GM-VFNS predictions for inclusive central and forward lepton production from the sum 
of charm and bottom hadron decay compared to ALICE ^-0] data. 

even over a largre pt range, up to 100 GeV. This agreement in this large pt range also shows 
that the evolution of the PDFs and FFs with increasing scale is very well accounted for by 
the data. 

The second group of data comes from the ALICE collaboration. These are cross sections 
da/dpr for muon production at forward rapidity in the region 2.5 < y < 4.0 at = 7 
TeV |7| and at \/S = 2.7 TeV The data for da/dpT are presented as a function of px 
beteen 2.0 < Pt ^ 12 GeV, integrated over the total rapidity range and separately over five 
bins: (1) 2.5 <y < 2.8, (2) 2.8 < y < 3.1, (3) 3.1 < y < 3.4, (4) 3.4 < ?/ < 3.7 and (5) 
3.7<y<4.0asa function of px in the same transverse momentum range . The data and 
the comparison with our predictions are shown for the total y-range in FigJH (second figure 
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= 2.76 and third figure y/S = 7 TeV) and for tlie cross sections integrated over tlie 
five y bins (1) to (5) in FigJH (most right figure) annd in FigJS] The agreement with the 
data is very satisfactory. The data points with their errors he mostly between the curve for 
the default scale choice (full line) and the prediction for the maximal scale choice (upper 
dashed line). In addition we show in FigJH also the results for the later published ALICE 
cross sections in the central rapidity region {\y\ < 0.5) (moost left figure). The agreement of 
the ALICE data j^l and our predictions with the default scale choice is again quite good. 
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FIG. 5: GM-VFNS predictions for inclusive forward lepton production from charm and bottom 
hadron decay compared to ALICE ?! data in various rapidity regions. 
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IV. SUMMARY 



We have calculated the cross sections for inclusive lepton production originating from 
heavy flavour decays at LHC cm. energies of 2.75 and 7 TeV in the framework of GM-VFN 
scheme and compared them with measured cross sections da/dpT in different rapidity regions 
of the ALICE, ATLAS and CMS collaborations at the LHC. We found good agreement with 
all the experimental data inside the experimental and theoretical accuracies for both, the 
lepton data coming only from bottom hadron decays (CMS and ALICE data) and the lepton 
data for the sum of charm and bottom hadron decays. This shows that within the considered 
kinematical as given by the LHC collaborations the description of the inclusive production 
of charmed hadron and bottom hadron is well accounted for. Except the CMS data all 
other inclusive lepton production data from the LHC collaborations have been compared 
to predictions of the FONLL approach 3| and similar good agreement between data and 
theoretical predictions has been found. 
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